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We give a brief overview of the way atomic physics is now developing in a way reminiscent of 
the optics revolution of the 1960's. Thanks in particular to recent developments in atomic trapping 
and cooling, the new field of atom optics is rapidly leading to exciting new developments such as 
nonlinear atom optics and quantum atom optics. We illustrate these developments with examples 
out of our own research. 
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I. INTRODUCTION 



Atom optics has witnessed considerable progress in the last few years. A number of optical elements, including 
i— h ■ atomic mirrors M and gratings H , have been demonstrated, and several typ es of atom interferometers have been 
built In addition, advances in laser cooling and trapping of atoms |t^ [ll| have led to spectacular developments, 

including the observation of quantized atomic motion in optical lattices ]l2| , p^ |, and, in the last few years, the 
demonstration of Bose-Einstein condensation (]Hj-fl6| and the realization of a primitive atom laser (lj],[l^] . Atom optics 
using condensates as sources has now become reality, as illustrated, e.g., by experiments involving the interference 
' of two condensates , the Kapitza-Dirac diffraction of a condensate off a standing- wave grating |ll| , or proposed 
, experiments in matter-wave phase conjugation | pp| and matter- wave amplification pl| , p2[ . 

The early experiments in atom optics considered low density samples, where atom-atom interactions are negligible 
ON , and the atoms in the beam behave independently. We call this regime linear atom optics. When atom-atom inter- 



jle becomes dependent on the presence of other 
the matter waves analog of nonlinear optics, 



■ actions become important, the dynamics of a given atom in the sam 
| atoms, and one reaches the regime of nonlinear atom optics p3|-^5| 
. which was pioneered by Peter Franken in 1962 p6| . 

In addition to making a distinction between linear and nonlinear atom optics, it is important to also separate 
the ray optics from the wave optics regime: In the first case, it is sufficient to describe the atoms as point particles 
following classical trajectories influenced by light fields. Much of laser cooling belongs to that category. In contrast 
the wave atom optics regime, which is normally — but not exclusively — associated with ultracold atomic samples, 
I ■ requires a proper quantum mechanical description of the center-of-mass motion of the atoms, including the effects of 
' O I matter-wave diffraction. 

Finally, a last distinction of importance is between what might be called the quantum and classical regimes of atom 
optics, again in analogy with the electromagnetic case: the quantum regime of both atom optics and optics is that 
regime where the effects of quantum statistics play a significant role. Similarly to optics, where both the classical and 
the quantum field are governed by Maxwell's equations, both the "classical" and the second-quantized Schrodinger 
, fields are governed by the Schrodinger equation. 

The quantum regime of matter- wave optics is conveniently described within the framework of second quantization 
5^ ! of the matter wave field, which is described "classically" by the wave function ip(r, t), so that ^(r, t) — ► ^(r, t). This is 
analogous to the optical case, where the electromagnetic field is quantized, E(r,t) — > E(r,t). It should be remarked, 
however, that in contrast to optics, where the quantization of the field introduces new physics, such is not the case 
when second-quantizing the Schrodinger field at the low energies considered here. In this regime, the total number 
of particles is of course conserved, since effects such as pair creation are ignored. Hence, for our purpose the second 
quantization procedure is merely a convenient book-keeping mechanism which automatically accounts for particles 
being added and removed from a given state with the proper quantum statistics. Note however that we intentionally 
distinguish the quantum from the nonlinear regime of atom optics, since reaching the first one requires that a condition 
on phase-space density be met, while the second one is related to density only. 

The goal of this paper is to illustrate the recent developments in quantum nonlinear atom optics with the help of 
three specific examples. We first recall in Sec. II how collisions are the main source of nonlincarities in atom optics, 
showing the parallels between this and the situation in nonlinear optics. Section III presents a matter-wave four-wave 
mixing process, the generation of a phase conjugate (time-reversed) matter wave in multicomponent condensates. 



* Dedicated to Peter Franken on his 70th Birthday 
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Section IV briefly reviews a proposed scheme of an atom laser, and Sec. V discusses a parametric process taking place 
in matter waves, the so-called collective atomic recoil laser. Finally, Sec. VI is a summary and outlook. 



II. THE ROLE OF COLLISIONS 



The analogy between nonlinear atom optics and nonlinear optics is quite profound: In conventional optics, effective 
nonlinear equations for the optical fields result from the elimination of the medium dynamics, while in atom optics, 
nonlinear matter-wave dynamics results from collisions, which are in turn an effective manybody interaction resulting 
from the (partial) elimination of the electromagnetic field. Indeed, at the most fundamental level of quantum electro- 
dynamics there is no such thing as two-body interactions between atoms; rather, the potentials that describe them 
are the result of a series of approximations whose validity depends on the precise situation at hand. 

One simple way to illustrate how this works is to consider the near-resonant dipole-dipole interaction inside a cavity. 
Consider two two-level atoms of Bohr transition frequency uj a , and located at positions Xi and X2, with |xi — X2I = x 
along the axis of a ring cavity of length L. This system is described by the Hamiltonian 

H = H a + H af +H f , (1) 

where the atomic Hamiltonian H a is the sum of the individual atomic Hamiltonians 



i=l,2 i=l,2 



Pi are the momenta of the atoms of mass M, with [x,-,Pj] = ihSij, and or and a_ are the atomic raising and 
lowering pseudo-spin operators for the i-th atom. The field Hamiltonian is 

H f=^2^ } n a i a ^ ( 3 ) 

where the composite index (i = {k„,^} labels both the wave number k„, with k n = uj n /c, and polarization i of the 
cavity mode of frequency u> n . The creation and annihilation operators a* and satisfy the usual Bose commutation 

relation [a^,ot/] = 8^1. Finally, the atom- field interaction Hamiltonian H a f is the sum of the individual atom- field 
electric dipole interaction Hamiltonians 

H?} = -di ■ E( Xi ), (4) 

where 

di = ed(a^ + a^) (5) 
is the atomic dipole operator, aligned along e and of magnitude d, and 

E(r) =z^£„[a M e lk "- r Q-ff.c.]. (6) 

n 

Thereby, £ n is the electric field per photon in mode /i, whose explicit form depends upon the choice of quantization 
scheme p7[p8| . For a running wave quantization scheme and periodic boundary conditions appropriate for a ring 
cavity we have in one dimension £ n = [hcj n /(2eoL)] 1 ^ 2 , where w„ = 2nnc/L. In the rotating-wave approximation, the 
atom-field Hamiltonian (Q) reduces then to 



II 



where 

g^iixi) = (£ n d/h)(e e ■ e) exp(ik„ • Xj). (8) 

We ignore in this section the kinetic energy part in the Hamiltonian (|2|), and consider the situation where one of 
the atoms is initially excited, the other is in its ground state, and the radiation field is in the vacuum state, 
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|^(0)) = \eg0). (9) 

In the rotating-wave approximation, the Hamiltonian (Q) conserves the number of excitations in the system, so that 
at time t its state vector can be expressed as 

m)) = hty)\eg0) + b 2 (t)\ge0) + £ b^t)\ggl^), (10) 

with 6i(0) = 1 and 6 2 (0) = 6 M (0) = 0. 

From the Schrodinger equation ih\?p(t)} — H\?p(t)}, the equations of motion for the various probability amplitudes 
involved are readily found to be 

^(*) = ^T J 9ntiv.ty), 

= -iA^ty) - glMt) - i&M*). (11) 

where i = 1,2 and A„ = oj„ — uj a is the detuning of mode fi from the atomic transition frequency. 

The simplest way to solve the set of equations ([ll]) is to neglect the effects of interatomic propagation. This 
is appropriate provided that x/c <C where T is the single-atom free-space spontaneous decay rate and x is 

an interatomic separation. Assuming that the atoms interact with a broadband vacuum and that the Born-Markov 
approximation holds, one then finds for the probability amplitudes involving excited-state atoms p9| 

b l , 2 {t)= l -[C + {t)±C„{t)], (12) 

where 

C±(i)=eKp|-i[n a ±ni 2 (a:)]t|. (13) 

Here we have introduced the single-atom and the two-body complex frequencies fi s and 12 (a;). The explicit form of 
Q s is 

Q 8 =T + iA s , (14) 

where 

t Jl roo j2 

T = 2Y / \g fli \ 2 6(c\k\-u a ) = — — — / dfcc|fc|7r*(c|fc|-w B ) = ^ (15) 
/* 

is a one-dimensional version of the free-space spontaneous emission rate, and 

A s = -r- / duj — — — (16) 

7T J U} a W - L0 a 

is the one-dimensional, two- level atom version of a Lamb shift. 

More interesting in the present context is the two-body complex frequency ^l\ 2 {x) is defined as 

n 12 (x) = r 2 (x) + iv dd (x)/h (17) 

and its real part T 2 (x) accounts for the modulation of the spontaneous decay rate of one of the atoms due to the 
presence of a second atom at a distance x, while its imaginary part is proportional to the dipole-dipole interaction 
potential Vdd{x). For the one-dimensional situation considered here, one finds explicitly 



n 12 (x) = I dkc\k\e lkx 



it5(c\k\ — u) a ) — iP 



1 



(18) 



so that 



T 2 (x)=T I duj—5{uj-uj a ) = Tcos(k a x) (19) 
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and 



Vdd {x) = -hT- r dco— = -hTsm(k a x), (20) 

K Jo U) a w-w a 

where k = tv/c and k a = tv a /c = 2Tr/\ a . For the more familiar three-dimensional case, the trigonometric functions 
appearing in these expressions are replaced by combinations of Bessel functions, and Vdd(r) falls off as 1/r for 
large interatomic separations, instead of being periodic [|30| -p2j. Note also that as a result of the Born-Markov 
approximation, these expressions neglect propagation and are independent of time. 

We see, then, that the two-body dipolc-dipole interaction between atoms results as advertised from the adiabatic 
elimination of the dynamics of the electromagnetic field. This is analogous to the nonlinear optics situation, except 
that the roles of the atoms and the electromagnetic field are reversed. 

In the following sections we discuss several examples of quantum nonlinear atom optics which show how the familiar 
ideas of nonlinear optics can be readily transposed to the new situation. We consider first four-wave mixing in a 
multicomponent Bose-Einstein condensate, which can lead, e.g., to matter- wave phase conjugation. We then discuss 
in Sec. IV a form of matter-wave "amplification" reminiscent of laser action in optics, and in Sec. V a parametric 
process, the low-temperature version of the collective atom recoil laser. Finally, Sec. VI is a conclusion and outlook. 



III. FOUR- WAVE MIXING: MATTER- WAVE PHASE CONJUGATION 



A particularly close analogy can be established between the dynamics of a spin-1 multicomponent condensate, as 
recently realized in sodium experiments and the situation of degenerate four-wave mixing in optics. Specifically, in the 
zero-temperature limit a g-component condensate can be thought of as a g-mode system, whereby the various modes 
are coupled by two-body (and possibly higher-order) collisions which result in the exchange of particles between these 
modes. 

In particular, for the case of a 23 Na condensate in an optical dipole trap one can achieve situations where the m = 
state is macroscopically populated while the m = ±1 states are weakly excited. One can then think of the first state 
as a "pump" or "central" mode, while m = ±1 form side modes which are coupled via the pump. This leads to 
familiar effects such as degenerate four-wave mixing and matter-wave phase conjugation. 

Consider then a condensate of 23 Na atoms in their F = 1 hyperfine ground state, with three internal atomic states 
\F = 1, m = —1), \F = l,m = 0) and \F = l,m = 1) of degenerate energies in the absence of magnetic fields. It is 
described by the three-component vector Schrodinger field 

*(r,t) = {*_!(!•,*), tf (r,t),*i(r,i)} (21) 

which satisfies the bosonic commutation relations 

[* i (r,t),*\(v',t)]=6 ij 6(v-r'). (22) 

Accounting for the possibility of two-body collisions, its dynamics is described by the second-quantized Hamiltonian 

H = J dr*t(r,t)J7 ¥(r,t) + J{dr}^(r u t)&(r 2 , t)V(n - r 2 )*(r 2 , t)¥(n, i), (23) 

where the single-particle Hamiltonian is 

H = p 2 /2M + U(r) (24) 

and the dipole trap potential U(r) does not depend on the hyperfine magnetic state m. 

The general form of the two-body interaction V{r\ — r 2 ) has been discussed in detail in Refs. p3]]34| |. Labeling the 
hyperfine states of the combined system of hyperfine spin F = Fi + F 2 by |/, to) with / = 0, 1, 2 and to = — /, 
it can be shown that in the shapeless approximation the two-body interaction is of the general form ]33] ] 

2 

V(r l -r 2 ) = 5{v 1 -r 2 )Y J hg f V fl (25) 
f=o 

where 

g f = 4mhaf/M, (26) 
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Vf = J2m l/; m )(/; TO l i s the projection operator which projects the pair of atoms into a total hyperfine / state and 
cif is the s-wave scattering length for the channel of total hyperfine spin /. For bosonic atoms only even / states 
contribute, so that 



V(n - r 2 ) = M(r a - r 2 )(g 2 V 2 + .90^0) = - r 2 ) (c + c 2 F 1 • F 2 ) . 



In this expression, 



c = 2(g + 2g 2 )/3, c 2 = 2{g 2 - g )/3. 
Substituting this form of V(ri — r 2 ) into the second-quantized Hamiltonian (p3|) leads to 



(27) 



(28) 



— + U(r) 
2M y ' 



* m (r,i) + ~ J dr{(c + c 2 ) [*+*t* x * x + *L 1 *L 1 *-i*- 



+2#J* (*t*i + *Li*-i)] + cq* Jtfjtf *o + 2(c - 02)^x^1^-1 + 202(^1*1^0*0 + H.c.)}. 



(29) 



This form of the Hamiltonian is quite familiar in quantum optics, where it describes four-wave mixing between a 
pump beam and two side-modes, which are identified with the field operators vE'o and *±i in the present situation. 

The three terms in the two-body Hamiltonian which are quartic in one of the field operators only, i.e. of the form 
^i^l^i^i, can be readily interpreted as self-defocussing terms, corresponding to the fact that the two-body potential 
is, for a positive scattering length and a scalar field, analogous to a defocussing cubic nonlinearity in optics. The terms 
involving two "modes", i.e. of the type ^t^^rt^-, conserve the individual mode populations of the modes and simply 
lead to phase shifts. Finally, the terms involving the central mode *o and both side-modes are the contributions of 
interest to us. They correspond to a redistribution of atoms between the "pump" mode *o and the side- modes *±i, 
e.g., by annihilating two atoms in the central mode and creating one atom each in the side- modes. This is the kind 
of interaction that leads to phase conjugation in quantum optics, except that in that case the modes in question are 
modes of the Maxwell field instead of the Schrodinger field. Note also that a similar mechanism is at the origin of 
amplification in the Collective Atom Recoil Laser (CARL) |35| |37| , see Sec. V. 

In the Hartree approximation, which is well justified for condensates at T = 0, the dynamics of a condensate 
described by the Hamiltonian (G9j) is governed by the system of coupled nonlinear Schrodinger equations 



i<j)-i(r,t) = - 



— + U(t) 
2M y ' 



0-i + NfafoMl + [(co + c 2 )(|<M 2 + l^ol 2 ) + (co - c 2 )|(/)i| 2 ] ( />-i} 



i</>o(r,i) = 
i0l(r,t) = 



P 

2M 

Pi 

2M 



U(r) 



U(r) 



<po + AT{co|0q| 2 0o + (co + Pa) (10- 



t> a + 2c20i0_i0o} 



0i + Af{c 2 0o0 o 0-i + [(co + c 2 )(|0i| 2 + |0o| 2 ) + (c - c 2 )|0_ x | 2 ]0i}. 



(30) 



Consider for example a situation where the central mode, described by the Hartree wave function tfio, is initially 
strongly populated while the side-modes <f>±i are weakly excited. In other words, we consider the phase conjugation 
of a weak atomic beam from a large condensate. It is then appropriate to introduce the matter-wave optics equivalent 
of the undepleted pump approximation, whereby 



0o-O, 

and the problem reduces to a set of coupled mode equations for the two side-modes 
a catalyst for the coupling between them. 

We take the trap potential U(r) to be of the harmonic form 

U(r) =MloI{x 2 +y 2 )/2 



(31) 

i±i, the central mode acting as 



(32) 



that is, we assume that the dipole trap confines the atoms in the transverse plane [x, y), but not in the longitudinal 
direction z. This geometry allows one to consider side-modes propagating along that axis, rather than bouncing back 
and forth in an elongated trap. In case of tight confinement in the transverse direction, we can assume to a good 
approximation that the transverse structure of the condensate is not significantly altered by many-body interactions 
and is determined as the ground-state solution of the transverse potential. 

Expressing the Hartree wave function associated with the hyperfine level m as 
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ro (r,i) = ip ± (x,y)(p m (z,t)e iuH,t , 



(33) 



we then have 



\-- 


( d 2 


d 2 \ 


2M 


\dx 2 ' 


h dy 2 ) 



(34) 



The physical situation we have in mind is that of a weak "probe" in the hyperfine state m = — 1 propagating 
toward a large condensate in state m = and at rest in the dipole trap, and generating a backward-propagating 



conjugate matter wave in the hyperfine state m = 
wave function as 



tp(z,t) 



(px{z,t) 



T. Hence we express the longitudinal component of the Hartree 

V>-i(z,i)e- ife2 
2-00 cos(kz) 
^i{z 7 t)e ikz 



where the slowly varying envelopes ip m of the Hartree wave function components m 
ities 



d 2 



< k 



9 I 



< k 2 \^ 



±1 satisfy the familiar inequal- 



(35) 



The "pump" wave function cpo is described by a standing wave, a configuration that can be achieved for instance by 
interfering two condensates fl8| ], in a grating matter- wave interferometer, or in the CARL j37|. To first order in 
the probe and signal fields, this geometry leads to a linearized system of two coupled- mode equations for the probe 
and condensate fields. In the stationary state and projecting out the transverse part of the wave function they reduce 
to 



hk d 

*2M = -^[2( c o + c 2 )poV>-i(z) + c 2 ^oV>i(-2)], 

hk d 



where po = \ipo\ 2 and 



J dxdy\(p ± (x,y)\ 4 
S dxdy\tp±(x,y)\< 



(36) 



(37) 



These equations are of course familiar from optical phase conjugation and their solution is well-known. The evolution 
of the phase conjugate wave ipi contains a term proportional to the density po of the condensate and the field itself. 
In the absence of the second term, it would simply lead to a phase shift of ipi- Physically, it results from the self- 
interaction of the conjugate field, catalyzed by the condensate (pump) component. Its origin can be traced back to the 
term proportional to \E , {^'{5'o^'o m the Hamiltonian (p9|). The second term, in contrast, couples the two side- modes 
via the condensate and is responsible for phase conjugation. 

The general solution of Eqs. (|36|) reads p8[ 



cos(|k|L) 



cos(|k|L) 



( - ie- lp Bin(|«|«)^(L) + cos(|k|(z - L))^-i(0)) 
( cob(\k\z)iI>i(L) + i e' lP sin(| K |( 2 - L))i>*_ x (0)) , 



(38) 



where a = 2Nrj(c + c 2 )p , k = S' , and e J/3 = k/\k\ 



hk/2M 

For the probe ^_i(0) incident at z - 
input plane z — becomes 



= and no incoming conjugate signal ipi(L) 
^(0) = -ie- i/3 tan(|«;|£)^* ^O), 



0, the conjugate wave in the 

(39) 



which demonstrates that the interaction of the probe and the condensate results in the generation of a counterprop- 
agating phase-conjugated signal. 

In addition to its interest from a nonlinear atom optics point of view, matter-wave phase conjugation could also be 
used as a diagnostic tool for Bose-Einstein condensates. For instance, we noted that the parameter \k\L is proportional 
to the difference in scattering lengths between the singlet and triplet states. Hence, this quantity could in principle 
be inferred from phase conjugation measurements. 
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IV. THE BINARY COLLISION ATOM LASER 



As a second example illustrating the close analogy between nonlinear optics and nonlinear atom optics, we now 
consider the atom laser. There has been quite a bit of confusion about what is meant by an "atom laser" in the 
past, apparently associated with the fact that since the number of atoms is conserved, one cannot possibly achieve 
their amplification.^] But of course, this is not the point. The main purpose of an atom laser is to create a coherent, 
macroscopic population in a given center-of-mass mode of atomic motion, such as to create a coherent atomic beam. 
Clearly, the number of particles in a given mode needs not be conserved: it can readily be amplified by cleverly 
extracting atoms from a large reservoir. 

A primitive atom laser has been demonstrated experimentally by the MIT group, which outcoupled a fraction of 
a sodium condensate from a magnetic trap via rf coupling between the weak-field-seeking mjp = -1 state and the 



strong-field-seeking mp = 1 state 17 lq|. We concentrate here on a different scheme, which has so far not been 



realized experimentally, the binary collision atom laser 40 - 44] . In this system, the matter- wave resonator consists 



e.g., of an optical dipole trap. In order to concentrate on the essential dynamics only three out of the multitude of 
trap modes are taken into account explicitly. 

The atom laser then operates as follows: Bosonic atoms in their ground electronic state are incoherently pumped 
into a trap level of "intermediary" energy (mode 1). There they undergo binary collisions which take one of the 
atoms involved to the tightly bound laser mode 0, whereas the other one is transferred to the heavily damped loss 
mode 2. This latter atom leaves the resonator quickly, thereby providing the irreversibility of the pumping process. 
A macroscopic population of the laser mode can build up as soon as the influx of atoms due to pumping compensates 
for the losses induced by the damping. 

In the description of this laser scheme one has to take into account that in addition to the pumping collisions 
other types of interatomic collisions can also occur. These considerations lead to an ansatz for the atom laser master 
equation of the form 

W = - = [H + H col , W] + K V[c }W + Kl {N+ l)V[ Cl ]W + k x NV[c\)W + k 2 V[c 2 ]W. (40) 

In this equation, we use the second quantized formalism in which each center-of-mass atomic mode is associated with 
a bosonic annihilation operator c;, and W denotes the atomic density operator .n The free Hamiltonian is given by 



H = hoJic\ci, 



u>i being the mode frequencies. The general form of the collision Hamiltonian can be written as 

Hcol = KVijkic\c^CkCi (41) 

i<j,k<l 

with Vijki the matrix elements of the two-body interaction Hamiltonian responsible for the collisions. It is sufficient 
for the present purpose to restrict our attention to the reduced form 

Hcol = M^02iiCo4 c i c i + ViW2c\c\c c 2 + V 00 o clclc c + V 0W iclc{c ci + Viinc{c\ciCi) (42) 

in which (besides the pumping collisions) only those collisions are retained which are expected to have the most 
significant influence on the phase dynamics. 

It is worth emphasizing that there is a fundamental difference in the way collisions are handled in conventional 
atomic physics and in the present situation. Usually, and for instance in the phase conjugation problem of the 
preceding section, collisions arc handled in terms of scattering amplitudes between initial and final states. In matter- 
wave resonators, however, such an approach becomes meaningless. Rather, the collisions are seen to induce transitions 
between cavity modes. Hence, the collision Hamiltonian ( [42| ) takes the form of a mode-coupling Hamiltonian between 
three cavity modes. It is therefore mathematically closely related to that of Eq. (|2^), except that the coupling is 
now between center-of-mass modes, rather than magnetic sublevels. Another difference is of course that in phase 



'Another problem is that the L in the acronym for laser stands for Light. But this is history repeating itself: in the early 
days of lasers, they were commonly called "optical masers." 
2 Since we consider ground state atoms only, they are fully described by their center-of-mass quantum numbers. 
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conjugation, we are interested in a coherent interaction, and neglect the effects of dissipation, which are known to be 
detrimental to the achievement of phase reversal. In contrast, the laser and the atom laser are open systems that rely 
explicitly on the presence of pump and probe mechanisms to achieve steady-state operation. 

The damping rates of the cavity modes are given by the coefficients Ki, and the strength of the external pumping of 
mode 1 is characterized by the parameter N, which is the mean number of atoms to which mode 1 would equilibrate 
in the absence of collisions. The superoperator T> is defined by 

V[a]P = aPa* - |(a+aP + Pa) a) (43) 

with arbitrary operators a and P. 

In order to achieve a sufficiently high degree of irreversibility it is necessary that K2 is much larger than the 
damping rates of the other modes. This suggests adiabatically eliminating this mode, an approximation that leads to 
the simplified master equation ^ ^] 

P = -jri H c P\ + KoV[c ]p + Ki(N + l)V[a]p + KiNV[c[]p + TV[clc\]p. (44) 

Equation ( fjif ) is written in the interaction picture with respect to Ho = Tiujqc^cq + hcdic\c\, and the reduced density 
matrix p is p = Tr m0{ j e 2 [W] . The reduced collision Hamiltonian H c is 

H c = ?i(VboooCo c o c o c o + Vmic[c|cici), (45) 

and r = 4|Vo2h| 2 /k. Consistently with Ref. jyj we call the limiting cases T <C ko and V ^ no the weak and strong 
collision regimes, respectively. The reduced master equation (^) forms the basis of most studies of binary collision 
atom lasers. 



The master equations (40) and (|44|) can easily be solved numerically using standard quantum Monte Carlo simulation 
techniques [ ^5|j46| . This is discussed in detail in Ref. jf3). Figure 1 shows the probability P(n) of having n atoms in 
the laser mode as a function of time for the case where only V0211 = V*i 02 is taken to be non-zero. The parameters 
used are ko = O.OIki, N — 2, k 2 = IOki and V0211 = 2ki, the time is in units of k^ 1 . We see that the population 
of the lasing mode builds up from zero to a state that suggests a Poissonian distribution with a mean number of 
atoms no = 40. In contrast, the evolution of the pump and decay modes (not shown) shows no significant build-up of 
population in steady state. 

Among the most important characteristics of a laser are the coherence properties of its output, in particular its 
linewidth. In order to obtain an analytical approximation for the laser linewidth in the two-mode system (Q) a 
linearized fluctuation analysis can be performed p8j , f47|] . Assuming as usual that the first-order correlation function 
Co( r ) — ( c o( r ) c o(0)) is determined only by the phase fluctuations one obtains in this way 

C (t) = rcoe"* e X p[-i^ (r)], (46) 
where the equilibrium populations of laser and pump mode are given by 

no = ~{N-n l ), (47) 

2 K 



fix = v/^Tr, (48) 

the threshold condition being N > \J kq/T. The time-dependent deterministic phase drift is indicated by 4>o- This 
phase drift leads to a shift of the center of the power spectrum (the Fourier transform of the correlation function) by 
an amount 2Vbooo^o + Vbioi^i with respect to the collisionless case. The behavior of the correlation function Cq(t) 
is thus essentially determined by the coefficient <j^ (r) which is the covariance of the phase fluctuations in the laser 
mode. 

This coefficient can be evaluated explicitly in a straightforward way, however the ensuing expression is rather 
complicated. In the following we restrict the discussion to two limiting cases which illustrate the essential aspects of 
the influence of the elastic collisions on the laser linewidth. (It should also be noted at this point that Co(t) does not 
depend on elastic collisions between pumping mode atoms, which are characterized by the parameter Vim.) 

(a) V0101 = 0. Equation ( f46| ) implies that the most important aspects of the correlation function can be inferred 
from the study of O0 O (r) in the time interval where it is smaller than or of the order of unity. Examining the behavior 
of C0 o (t) as a function of Vbooo (with all other parameters kept constant) one can distinguish between two different 
regimes. For small Vqooo the time evolution of C0 o (t) relevant for Cq(t) is well approximated by 



8 



a^ir) ~T[w + K /(2n )}, (49) 

where 

™ = Voooo^? 2 + 2VWT+-^= ) . (50) 



1 o 



at- vWr 



For such values of Vboooi Co( T ) decays therefore exponentially, and the power spectrum is Lorcntzian. If w Ko/(2no) 
the linewidth is proportional to VQ OOQ ho. In contrast to the situation with conventional lasers, it increases linearly 
with the number of atom s in the laser mode. 

In case Vbooo > J ^ "1^ w ^ n w ' = W^oooo we nn d a different behavior of the correlation function. Under these 
circumstances it decays like a Gaussian, i.e., 

ct 0o (t) ~ W* 000 n T 2 . (51) 

The spectrum is thus itself also of Gaussian shape and its linewidth proportional to Vbooo \fno- The atom laser 
linewidth still increases with no, albeit less dramatically than in the preceding case. 
(b) Vbooo = 0. In this case the expansion of <T0 o (t) to leading order in uq yields 

<Wr) = (Mgi + £-) r (52) 
\ 2rn 2n / 

The correlation function thus decays exponentially for all values of Vbioi- A qualitative change in behavior as in 
the previous situation does not occur. The linewidth of the spectrum is now proportional to Vo 2 101 /no, and becomes 
narrower when the population of the laser mode is increased, v ery much like the familiar Shawlow-Townes linewidth 



of conventional lasers. The analytical estimates of Eqs. (|49j), (51), and ( [52] ) are in good agreement with numerical 
quantum Monte Carlo simulations E3]. 



V. PARAMETRIC AMPLIFICATION: THE COLLECTIVE ATOM RECOIL LASER 

As a final example of nonlinear atom optics in ultracold atomic systems, we briefly review some of the most salient 
aspects of the ultracold atoms operation of the Collective Atomic Recoil Laser, or CARL |35|. This device consists 
of three main components: (1) the active medium, which consists of a gas of two-level atoms, (2) a strong pump 
laser which drives the two-level atomic transition, and (3) a ring cavity which supports an electromagnetic mode (the 
probe) counterpropagating with respect to the pump. What makes the CARL interesting is that the initial state, 
consisting of a thermal cloud of atoms and no photons in the cavity, is exponentially unstable. Laser oscillations appear 
spontaneously in the probe mode correlated with the appearance of a density modulation in the atomic sample. The 
original purpose of the CARL, which operates then at room or higher temperatures, was the generation of a tunable 
coherent light field from atoms in a way similar to light amplification in free-electron lasers, i.e., gain correlated with 
bunching and in the absence of population inversion. But as we shall see, when operating at ultracold temperatures 
it can simultaneously parametrically amplify spatial side-modes of a Bose-Einstein condensate. This is the aspect of 
the CARL that we concentrate on in this section. 

In the absence of collisions, the second-quantized Hamiltonian of a sample of two-level atoms interacting with a 
classical pump laser and a counterpropogating probe cavity mode is 



H = Y^H(k) + hcqA^A. 



(53) 



Here H (k) is given by 



H(k) 



2 1,2 



■c\{k)c g {k) 



2 1,2 



2m 



hci! I 4(fc)c e (fc) 



2m 

%-±e iULt c\{k + k L )c e {k) + ihgA^cUk - q)c e (k) + H.c. 



(54) 



where the bosonic matter wave operator c g (k) annihilates a ground state atom of momentum Tik, and c e (k) annihilates 
an excited atom of momentum hk, with 
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l (k),cUk')} = [c e (k),cl(k , )}=5 kk ,, (55) 



all other matter- wave commutators being equal to zero. The pump is characterized by its frequency cjl, its Rabi 
frequency O^, and its wavenumber k^. The probe field operator A annihilates a photon with wavenumber q, and the 
atom-probe coupling constant is g = d[cq/ (2eoLS)] 1 ^ 2 , where d is the atomic dipole moment, and LS is the cavity 
volume. Note that the atomic recoil is explicitly included in the electric dipole interaction coupling the electronic 
ground and excited states. 

The Hamiltonian (^ij) readily yields the Heisenberg equation of motion for the probe field mode as 

^A = -icqA + gJ2 c U k -l)c e (k). (56) 

k 

Hence, all that is required to determine the field evolution arc bilinear combinations of atomic creation and annihilation 
operators. Their equations of motion are in turn 

j t cl(k)c e (k') = l[H,cl(k)c e (k% (57) 

Because the Hamiltonian (|34|) is bilinear in the atomic operators, we observe readily that the equations of motion 
([57]) involve only bilinear combinations. This is a direct consequence of the fact that collisions are neglected in this 
model. In the language of manybody theory, this means that the BBGKY hierarchy is exactly truncated. If the probe 
field is treated classically, then Eq. (|5^) becomes a c-number equation involving only the matter-wave expectation 
values (c g (k — q)c e (k)), which are nothing but the single-particle atomic density matrix elements. Since these elements 
are coupled only to other single-particle density matrix elements, the system can be solved without any truncation 
scheme being required, albeit numerically in general. If instead of treating probe field classically we had chosen to 
treat the atomic single-particle density matrix classically, we would have reached the same result. Thus converting 
Eqs. (|5^) and (|57|) into c-number equations is equivalent to a classical field theory for both the atomic and optical 
fields. This approach is discussed in detail in Ref. J37[ , and we do not pursue it further here. 

Since we are interested in the interaction between light and a Bose condensate, we must be careful that electromag- 
netic heating does not occur. Hence, the optical fields must be far off-resonant from any electronic transition, and 
the upper electronic levels can be adiabatically eliminated. The atoms are then described as a scalar field, since only 
their electronic ground state remains. In that case one finds the system is described by the effective Hamiltonian 

H= ^Y, fc2 4( fc K( fc ) + hc ^ A + £ [^Le-^A^cKk - 2k Q )c g (k) - H.c] 



A L \ 4 



g\ 2 A^A )^2cl(k)c g (k), (58) 



where we have introduced the recoil kick 2fco = q + ki,, and the pump detuning = cj^ — cj q . 

If the atomic sample is initially a condensate at zero temperature, the dominant mode, at least for short times, is 
the k = mode. It is macroscopically populated, while all other matter wave modes are in a vacuum. The equation 
of motion for the condensate mode is 

^c s (0) = ^- [Clle-^Alcgpko) - fi L e^^ Cg (-2fc )] (59) 

This equation indicates that the condensate mode is coupled via the electromagnetic field to the two side-modes at 
±2fco. In addition, each side mode is coupled to its neighbors at intervals of ±2/co- This is mathematically similar to 
the situation of Sec. Ill, see Eq. (|30|), with, however, three major differences: (a) in the case of phase conjugation, 
the mode coupling is between different magnetic sublevels of the atomic field, while it is now between momentum 
states; (b) the coupling is now due to a Raman-like coupling induced by two counterpropagating fields, instead of 
ground-state collisions. Hence the mode coupling now involves an infinite manifold of matter-wave modes and two 
electromagnetic fields modes — one of them treated classically in the present example — instead of four matter-wave 
modes. This is of course a direct consequence of the fact that collisions are the result of a phenomenological approach, 
eliminating the electromagnetic vacuum as we have seen in Sec. II. (c) Instead of three magnetic sublevels, we now 
have an infinite hierarchy of center-of-mass modes of the matter-wave field to deal with. 

This system is similarly related to the atom laser of Sec. IV, except that in that case, matter-wave modes are 
selected by the atomic resonator, as well as possibly by the selection rules of the collision Hamiltonian M, whereas 
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here they are selected from the continuum of modes simply by conservation of momentum. In addition, there are no 
atomic pump and decay mechanisms in the present model, which describes a fully coherent interaction. 

Let us now discuss this hierarchy in some more detail. We note that each mode k is directly coupled only to its 
neighboring modes k ± 2ko- But except for the condensate mode k = 0, all modes are initially empty, so in the early 
stages the dominant dynamics results from the coupling between the condensate mode and its two neighboring modes. 
Neglecting then the higher-order modes, and further treating the condensate mode as a constant c-number, a sort of 
undepleted pump approximation for a classical atom-laser field and an excellent approximation at T — and for a 
sufficiently large condensate, we find that we have reduced the system to a linear three-mode problem. It is easily 
shown that this reduced problem can be described by the effective Hamiltonian 

H = Ahu>R fclc_ + c^c + — Aa'a + \ o^cl + a) c + + c^_a + c^a ^ , (60) 

where ujr = hk^j2m, is the atomic recoil frequency, 

c ± = e 4|Oi|2 ' /4Ai c s (±2fc ), (61) 

a = -tG^Ai/lflliniHAiDe^A, (62) 
x = MM^, (63) 

A=^, (64) 

u) = cq— \g\ 2 N/ Al and N is the mean number of atoms in the condensate. We see that x 2 is an intensity parameter, 
proportional to the product of the intensities of the pump laser and the initial condensate, and A is simply the 
pump-probe detuning in units of 4lJr. The Hamiltonian ( |60| ) gives the full quantum field theory description of the 
zero temperature CARL, and is valid for all times short enough so that (c±c±) -C N and (a^a) <C |£1l| 2 /4M 2 - 

The presence of terms such as atc\_ in Eq. (60) immediately brings to mind the non-degenerate optical parametric 




amplifier |52j] which generates highly non-classical optical fields. These fields exhibit two-mode intensity correlations 
and squeezing, and have been extensively employed in the creation of EPR photon pairs for fundamental studies of 
the laws of quantum mechanics. Our system represents a generalization in that we now have three entangled quantum 
fields, and is especially interesting in that two of the fields are atomic rather than optical. A detailed analysis of 
quantum mode coupling in the CARL can be found in Ref. pl| . Here we focus on the intensity dynamics and 
fluctuations. 

The dynamics of the system can be determined by solving the three coupled mode equations 

(65) 

\c + (t)J 

where we have introduced the dimensionless time variable r = AujRt. 

The spectrum of (|6^) has been studied in detail in |H]|57]], with the result that under certain threshold conditions 
the eigenvalues take the form Ai = u>%, X 2 = £1 + iT, and A3 = Q — iT, where u>%, f2, and T are all real quantities. 
This means that the time-dependance of the operators will, after some transients, grow exponentially in time at the 
CARL growth rate T. 

If we now consider an initial state where the atomic side modes are in the vacuum state and the probe mode has 
been injected with an initial field in a coherent state |a), we find that after the transients have died away the mean 
intensities of the three modes are given by 

I A (r) ^ (a^r)a(r)) » (\a\ 2 R 2 n + R 2 l2 )e 2T \ (66) 
J_(r) = (cL(t)c_(t)) « (\a\ 2 R 2 1+ R 2 2 )e 2r \ (67) 

and 
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J+(r) = (cl(r)c + (r)) « (|a| 2 i& + i^Je 2 ^. (68) 

Here the coefficients {Rij} are given by Rij — Iv&Vqj |, where the columns of the matrix u are the eigenvectors of the 
3x3 linear system described by Eq. ( |65| ) . For a given set of control parameters x an d A these coefficients are simply 
constants, whose analytic expressions, while straightforward to derive, are too unwieldy to reproduce here. Thus we 
see that the intensities have a stimulated component, proportional to |a| 2 , and a spontaneous component, which is 
present even when all three modes begin in the vacuum state. In this case the system is triggered by quantum noise 
in the form of fluctuations in the atomic bunching. 

The second-order equal-time intensity correlation function is defined, e.g. in the case of the probe field, as 

[2] , v _ (aHT)a^T)a(T)a(T)) 
9 () ~ (at(r)a(r)) 2 ' [ > 

This gives a measure of the intensity fluctuations, and hence the coherence properties of the various modes. After the 
transient regime, the correlation functions for all three modes are given by 

^"'- g^biF (70) 

where the fluctuation function f(x, A) = \v^ 2 /« 31 | 2 is approximately unity near the region on the A axis where the 
exponential growth rate T is maximized for fixed x, and steadily increases as one moves away from this region. We see 
that in the spontaneous case (|a| 2 = 0) the intensity correlation functions are equal to 2, which is characteristic of a 
thermal, or chaotic field. However, as \a\ 2 is increased, approaches 1, which signifies a coherent, or Poissonian field. 
Thus we see that by varying readily adjustable experimental parameters, such as the injected probe field strength and 
frequency, the condensate atom number, and/or the pump intensity and detuning, we have the capability to vary the 
intensity fluctuations of the generated fields continuously between thermal and coherent limits. 



VI. OUTLOOK 



In summary, then, we see that atom optics is progressing along a path that closely parallels that of optics following 
the invention of the laser. In many situations, it is possible to understand an atom optical effect by simply reversing 
the roles of light and matter, and indeed, much of the inspiration leading to atom lasers and nonlinear atom optics 
results from such an approach. The pioneering work of Peter Franken and his colleagues play a central role in these 
developments. But what would probably intrigue Peter most is not so much the way in which atom optics and optics 
are similar, but rather those important aspects in which they differ, as results in particular from the fact that atoms 
have an internal structure, are massive, and are composite particles. Work along these lines will no doubt keep him 
amused for many years. 
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FIG. 1. Probability P(no) for having no atoms in mode as a function of time. The parameters for this plot are ko = .01, 
«i = 1, N = 2, k2 = 10, V0211 = 2, and Vbooo = Vbioi = Vim = 0, all given in units of k^ 1 . 
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